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Decaying and periodically kicked turbulence are analyzed within the GOY 
shell model, to allow for sufficiently large scaling regimes. Energy is trans- 
fered towards the small scales in intermittent bursts. Nevertheless, mean field 
arguments are sufficient to account for the ensemble averaged energy decay 
, E(t) ~ t~ 2 or the parameter dependences for the ensemble averaged total 

CN \ energy in the kicked case. Within numerical precision, the inertial subrange 

O ' intermittency remains the same, whether the system is forced or decaying. 

Q 

\D ■ I. INTRODUCTION 

Q | Three dimensional turbulent flows are characterized by a highly chaotic and intermittent 

transfer of energy from the outer length scale down to the dissipative, inner length scale. 
In this paper we will focus on two different kinds of turbulent flows: Decaying turbulence 
and periodically kicked turbulence. Both types of flow have already been analysed within a 
mean field theory ]T[EJ, but here we would like to focus also on intermittency effects which 
cannot be described in the mean field approach ||. 
> ; To be more specific: With decaying turbulence we mean a homogeneous and isotropic 

turbulent flow for which forcing is ceased from some time to on. Therefore, eventually all 
the energy will be damped because of dissipative effects, but the statistical properties of the 
decaying turbulent field are a priori not clear. With periodically kicked turbulence we mean 
q ! a turbulent field forced with short and very strong periodic pulses. 

The motivation of the paper is to study the effect of non-trivial forcing on the properties 
of the turbulence. Most numerical or theoretical studies assume a Gaussian random noise 
forcing, acting on the largest scales only. But for most practical flows the forcing protocol 
is obviously more complex and often periodic, be it pulsed flow through a pipeline or the 
earth's atmosphere driven by the periodical heating through the sun. Another example of 
periodically kicked turbulence is the numerical realization of homogeneous shear flow @] 
where periodical remeshing is necessary. In between the kicks the turbulence is supposed to 
be freely decaying. Therefore an understanding of decaying turbulence is required. Beyond 
this, decaying turbulence is of course one of the classical examples of turbulent flow and an 
extended literature exists, see e.g. |5|-^§. 

The aim of the paper is to explore the statistical properties of turbulence both in the 
decaying and kicked case. Therefore, we must have excellent statistics. This turns out to be 
prohibative in a direct numerical simulation (DNS) of the Navier-Stokes equation and that 
is why we revert to study the problem in the context of the GOY shell model of turbulence 
|10|| in which a scaling regime of many decades can be achieved. To our knowledge the 



> 



present study is the first on the classical problem of decaying turbulence with the help of a 
shell model. 

Shell models are defined by a set of hierachically coupled ODEs for the velocity modes 
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which try to reproduce the physics of the energy flux of Navier-Stokes equation. The dy- 
namical equations for the GOY models read [10-13 



u n+ iM n _i + c„n n _in n _ 2 )* + g5 n ,o (1) 

were n = 0, ... ,N — 1, k n = 2 n ko, a n = 1, b n = —5/2 and c n = —(1 — <5)/4 and the 
boundary conditions are a^-i = ajv = fejv = &i = ci = c 2 = in order for the GOY model 
to conserve the energy in the unforced and inviscid case (g — v — 0). Traditionally, the free 
parameter 5 is chosen to be 5 = 1/2. The values of the other parameter used were N = 22, 
g = (1 + i) • 10" 2 , and v = 1(T 6 . 



II. DECAYING TURBULENCE: GLOBAL PROPERTIES 

The stationary (i.e. forced) simulations of eq. (|T|) were performed using fourth order 
Runge-Kutta with viscosity explicitely integrated. For the decay run we used the same 
algorithm but increased the time step keeping it 1/10 of the dissipative time scale. During 
the decay process indeed the dynamics becomes slower and slower and investigation of long 
time properties is only possible using a scheme with an adaptive time step. 

We performed a very long stationary simulation of equation (|l|) from which we collected 
an ensemble of starting configurations (~ 2500 independent runs) for studying the free decay 
of GOY turbulence. The starting configurations of the ensemble were collected after some 
eddy turnover time in order to be statistically independent. We let each of the starting 
configurations decay according to equation ([!]) (i.e. with g = and all other parameters 
unchanged) . 

During the decay stage we measure the total energy E(t) = J2n=o \ u n{t)\ 2 as a function 
of the decay time t (the time elapsed from when we switched off the forcing). In figure |]the 
decay of the total energy on a single run is shown, i.e., no ensemble average is performed here. 
The total energy decays in burst, i.e., it is constant for some time and than very suddenly 
drops. As this behavior occurs on all scales, the step-like structure looks self similar. The 
averaged (over the whole decay) decay exponent is close to —2. For very large times it 
decays exponentially. The derivative of the energy decay is shown in figure || in order to 
better highlight the bursty structure. Also from figure ^ one can immediately realize that 
the typical time scales during the decay becomes larger and larger. The bursty structure 
of the energy decay in the GOY model has been analyzed in detail by Okkels and Jensen 



14],|15| and it reflects the intermittent behavior of the energy flux downscale. 

We now ensemble average E(t) by collecting various starting configurations as described 
above and letting them decay. In figure |3] we show the time decay of the ensemble-averaged 
total energy. The step like behavior of a single realization of E(t) is now completely smeared 
out and the decay law (E(t)) ~ t~ 2 is revealed. For large times the averaged decay is of 
course again exponential. 

We now set up a simple model which is able to describe the average energy behavior 
during the decay process. The model closely follows the mean-field model of reference [|TJ. 
The major assumption we make is to suppose that all the energy is contained in the smallest 
shells corresponding to the largest scales. For simplicity we assume that it is only in the 
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zeroth shell of eq. (J!]), E ~ (|wo| )• During the first part of the decay energy will disappear 
from the zeroth shell by being transferred to smaller scales 
£(E{t)) 3/2 with the solution 



Uq\ 



and hence 



(E(t)) 



-1/2 



(2) 



Notice that asymptotically this means that, for t ^> to, we expect the energy to decay 
quadratically in time, (E(t)) ~ t~ 2 , just as seen in figure [3| The reason that the mean field 
argument works for the averaged decay is that the required time for the transport of pulses 
downscale is determined by the large scale dynamics, see e.g. figure lb of Schorghofer [TJ[ . 

As energy is removed from the system the effective Reynolds number will decrease further 
and further. In particular there will come a time, t\, for which even on the zeroth shell the 
dissipative term will dominate with respect to the non linear term in equation ([J). From 

that moment on, the equation for the energy decay will be 



Wo 



-vkl 



and 



hence 4 (E(t)) ~ — vk\ (E(t)) whose solution is an exponential damping, 



(E(t)) = E(ti) ■ exp -2ukl (t - ti 



(3) 



As seen from figure ^, equations (0) and 
term behavior of the energy decay. 



correctly describe the respective short and long 



III. INTERMITTENCY IN DECAYING TURBULENCE 



The time behavior of an individual decay process is complicated because of the bursty 
structure of the decay, see figure ||. The presence of bursts is an essential feature of inter- 
mittency. During the decay process both the intensities of the bursts decrease and their 
duration increase orders of magnitudes. 

The question now is if or not the statistical properties of the turbulent fluctuations 
remain the same during the decay process. To answer this question we study the k-scaling 
of higher order moments. Rather than calculating the scaling of moments of the velocity 
u n itself, which within the GOY model shows unphysical period 2 and period 3 oscillations 
||13|| , we focus on the scaling of energy flux moments [13] E 



En,p(^) 



J (u n U n+ iU n+2 + ~U 



(4) 



Here X denotes the imaginary part. The E„ iP are free of the period 2 and period 3 artifacts 
in the spectrum and show very clean scaling properties ][H|]. With the angular brackets (. . .) 
we again denote ensemble average conditioned to a given decay time t. However, to obtain 
better statistics, here we had in addition to average over a short period of time. We chose 
a tenth of a decade and the given times refer to the end of that small time interval. 

The numerical results are shown in figure |] where we have plotted various moments of 
the fluxes as a function of decay time. One evident feature is the decrease of the Reynolds 
number as the decay goes on. This results from the strong increase of the dissipative scale 
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and consequently the shortening of the inertial range. In the remaining inertial subrange 
(ISR) the slope looks very similar, but whether it really is the same cannot yet be judged 
from this type of plot. 

Therefore, in order to explore whether the scaling in the forced and in the unforced cases 
are really the same, we calculate the ratio between powers of the fluxes computed at different 
decay times, see figure |5|. Though in the forced case the forcing is limited to the zeroth shell, 
slight deviations spread over the first three shells or so. However, in the ISR the scaling 
properties and thus the intermittency really seem to be the same, within statistical errors. 
Note that getting this conclusion from DNS would be very hard, due to the small extension 
of the inertial subrange. 

This result on the lack of dependence of intermittency on the forcing resembles analogous 
conclusions from analysing the effect of the viscous subrange on the ISR scaling: For both 
the GOY model and also for DNS She and coworkers 07 01 found the same degree of 
ISR intermittency, independently of whether normal or hyperviscosity was employed. Only 
slightly beyond the onset of turbulence there may be a small dependence on the type of 



viscosity within the GOY model [19 



IV. KICKED TURBULENCE 

Kicked turbulence has been analysed within the framework of a mean- field theory ||. 
In the kicked case the GOY model is forced with a delta-like (in time) forcing gd ni o with 
frequency /. With delta-like forcing we mean that the forcing g5 n fl is periodically turned 
on for a small time Atk ic k <C The presence of the forcing will sustain the energy flux. 
The turbulence level achieved in kicked turbulence depends on both the forcing strength 
A = gAtkick and the forcing frequency /. We employ the GOY model dynamics to explore 
this dependency beyond the simple mean field approach of reference [[|. 

The qualitative behavior of the energy, as a function of kicking frequency can be seen 
from figures |]. In figure |]a we can see the energy behavior in the laminar regime. After each 
kick the system jumps at (almost) constant upper levels. From this value it then decays for 
a long period before another kick is applied. 

In figure ||b we are in the transition regime towards turbulence. The time between 
two kicks has been decreased as compared to Fig. ||a and hence the system does not have 
sufficient time to fully relax. The kicks still heavily influence the macroscopic behavior. 

In figure |6]c the kicking frequency has been further increased and we are close to the 
turbulent regime. We see regions where the energy starts to pile up for a while, after which 
it relaxes very fast through bursts of energy. The macroscopic energy behavior is described 
by a competition between the kicking (with timescale 1//) and the energy decay between 
the kicks (with the time scale of the large eddy turnover time). 

In figure |6]d we are in the fully turbulent regime. Here the individual kicks are no longer 
important for the macroscopic energy behavior and they acts as a sort of average constant 
forcing. Energy can build up over many kicking timescales before it is released through an 



energy burst, whenever the phase relation is appropriate [14|. This behavior is the analog to 



the bursty structure of figures [TJ and in the decaying case: Energy can only be transported 



downscale if the phase relations happen to be appropriate [14]]. If not, an energy plateau 



forms in the decaying case or the energy piles up in the kicked case. 
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To get the dependence of the total energy just before the kick (Ei) and just after the kick 
(E u ) as a function of the kicking frequency and the kicking strength we again performed 
ensemble averages of many realizations. In figure [7] we show Ei and E u as a function of / for 
three different forcing strengths A. Just as in the mean field case || three regimes can be 
seen: (i) A laminar regime in which the upper energy level is constant and corresponds to the 
total energy input during the kick and where the time between the kicks is sufficiently large 
so that the energy completely decays; (ii) a transitional regime; and (iii) a turbulent regime 
where the average upper and lower energy level are equal as the forcing is experienced as a 
continuous forcing. The energy levels roughly scale with f l . In the last regime the features 
are different from those of the third regime in the mean field theory where the average 
energy of the lower level is always less than that of the upper level. The reason is that in 
the dynamical model energy can build up over many forcing periods due to phase blocking. 
Therefore, it is experienced as a continuous forcing. 

Nevertheless, at least the scaling can be obtained by a similar argument as employed 
in the mean field theory ||: If the energy of the GOY can be (roughly) approximated by 
the energy contained in the largest shell, E\ ~ (|uo| 2 ), then directly after the kick it will be 
E u ~ (| Mo + 9^tkick\ 2 ) = (| Wo + A\ 2 ) = Ei + 1A*jE~i + A 2 . Between two kicks the system is 
freely decaying. Therefore we can apply in between two kicks the mean field result of refs. 
]l],[2|] to connect E u and E x . If the decay starts from an energy value E u , after a time 1/ / , 
when the next kick is applied, the energy will be decayed to E[. The two energy levels are 
connected through the equation |],|2| 

-!- = 3 [F{Re{E{)) - F(Re(E u ))] . (5) 

JT 

Here the function F(Re) is defined as |l 

' 7 + Vi 2 + Re 2 ' 



Re 



(6) 



with 7 = 9. The Reynolds number and the energy in the GOY are connected by Re(E(t)) = 

\j\ L ^v^ ~- Here, L is the integral scale and r = L 2 jv a viscous time scale. 

Solving equation (H) we can find the value of the upper, E u , and lower, Ei, energy levels 
as a function of the forcing strength A and frequency /. In figure |8] this dependence is 
plotted, revealing the basic features as in the numerical figure [7[ For large Re 3> 7, we have 
F(Re) — 1/Re and one finds E u ~ E\ ~ /, just as in the numerical case. 



V. CONCLUSIONS 

We summarize our main findings: 

The GOY model is employed to study both decaying and periodically kicked turbulence. 
Energy is transfered towards the small scales in intermittent bursts, leading to stepwise 
behavior in the decaying case or energy pileups and subsequent bursts in the kicked case. 
In spite of this intermittent behavior, mean field arguments as developed in refs.[|I],[| are 
sufficient to account for the ensemble averaged energy decay (E(t)) ~ t~ 2 or the parameter 
dependence E u ~ E\ ~ / for the ensemble averaged total energy in the kicked case. The 
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reason that mean field arguments work here is that the decay and the kicking is determined 
by the large scale dynamics. 

For what concerns the statistical properties of decaying turbulence, our finding supports 
the idea that decaying turbulence has the same intermittency as stationary turbulence. In 
particular this finding allows to conclude that at least in the GOY model ISR intermittency is 
independent on the forcing mechanism of turbulence. The only relevant dynamics during the 
decay process seems to be a shortening of the inertial range (decrease of Reynolds number) 
but leaving the same intermittency properties. 
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FIG. 1. Time decay of the total energy E(t) as a function of time t for one particular realization. 
Note the step-wise, self-similar behavior. The average slope —2 is also shown. 
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FIG. 3. Ensemble averaged energy as a function of decay time. In order to give an idea about 
the statistical error, four sets are shown (shifted in y-direction). Each one is obtained by averaging 
500 uncorrelated decay processes. 
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FIG. 4. Second order flux moment £ Hj 2 at six different decay times. The upper curve is the 
original situation, the following curves reflect the situation one, two, three, four, and five decades 
in time later. Besides ensemble averaging over 500 ensembles, the respective data are also averaged 
over one tenth of a decade in time. 
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FIG. 5. Ratio of T, n> Q(t) to £ n ,6(io) one, two, three, four, and five decades after to when the 
decay started. For clarity the curves are again shifted in y-direction. 
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FIG. 6. (Top-Left): Total energy in the GOY as a function of time, in the laminar regime 
(/ = 6 • 1(T 6 , A = 0.01). (Top-Right): Closely after the laminar regime (/ = 1 • 10" 4 , A = 0.01). 
(Bottom-Left): Close to the turbulent regime. (/ = 1 • 10~ 3 , A = 0.01). (Bottom-Right): In the 
turbulent regime. (/ = 1.5 • lO^ 1 , A = 0.01). 
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FIG. 7. Numerical computation of the energy levels. E\ (lower curves) and E u (upper curves) 
as a function of the frequency / for three different kicking strengths A. For large / the slope is 
roughly E u ~ E t ~ / '. 
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FIG. 8. Energy levels E u (upper curves) and E\ (lower curves) as a function of the frequency for 
three different kicking strength, as they follow from the dimensional argument given in the text. 
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